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We study the competition between the spin-Peierls and the antiferromagnetic ordering in disordered 
quasi-one-dimensional spin systems. We obtain the temperature vs disorder-strength phase diagram, 
which qualitatively agrees with recent experiments on doped CuGeOa. 
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PACS numbers: 64.70Kb, 75.10Jm 

The discovery of the first inorganic spin-Peierls (SP) 
material CuGeOa opened the possibility to study the in- 
fluence of doping on the SP transition |l[ . The recently 
obtained phase diagram of doped CuGeOa has several 
surprising features |y-|5| . It turns out that doping, while 
suppressing the SP state, at the same time induces long- 
range antiferromagnetic (AF) order, with the Neel tem- 
perature initially increasing with the doping concentra- 
tion. Furthermore, a doping range is found where SP and 
AF order coexist. 

At first glance, it seems very strange that disorder 
(doping) may lead to the enhancement of some order pa- 
rameter (in this case the AF one). Also, the coexistence 
of the dimerized SP state, in which spins are bound into 
singlets, with a spontaneous subblattice magnetization 
that requires the presence of free spins, is rather puz- 
zling. In this Letter we address both these issues and 
obtain a phase diagram that is very similar to the ex- 
perimental one. We also predict a new phenomenon: a 
re-entrance transition from the dimerized SP state back 
into the undimerized state with decreasing temperature. 

Theoretically, the possibility of long-range magnetic 
order in doped SP systems was discussed in Refs. [0-^, 
where the lattice was treated classically and it was as- 
sumed that impurities "cut" the spin chains into finite 
segments. It was argued that the lattice relaxation in 
these segments results in the appearance of regions with a 
suppressed dimerization (close to impurities in the model 
of Ref. 10], or centered at kinks in the lattice dimeriza- 
tion in Refs. |§,^). The AF correlations that develop in 
these regions may, in principle, stabilize an inhomoge- 
neous state in which the SP and AF orders coexist. The 
enhancement of the magnetic susceptibility by disorder- 
induced kinks was also discussed in Ref. |1^]. Although 
these considerations provide a qualitative understanding 
of the magnetic ordering in doped SP materials, the de- 
scription of the thermodynamics of the mixed SP-t-AF 
state within the same approach is complicated and so far 
has not been given. 

In this Letter we consider a model that does allow for a 
detailed study of the competition between the SP and AF 
phases in the presence of disorder. Instead of considering 
disorder that randomly cuts chains into finite segments. 



we assume that doping results in small fluctuations of the 
spin-exchange constants on many bonds. Furthermore, 
we treat the lattice and the interchain interactions in 
the "chain mean fleld" approximation |11[] . Then, the 
effective single-chain Hamiltonian reads: 

Hs = 2_^ '^n.n+l^n^n+l — '^ /_^(^)"'S'^ i (1) 

n n 

where the exchange constants have the form: 



Jn,n+1 = Jo + (-)"A + SJn. 



n+1- 



(2) 



Here, A is the average value of the SP dimerization, and 
SJn,n+i is the random contribution due to doping. The 
antiferromagnetic order parameter h in Eq. (M is the am- 
plitude of the alternating magnetic fleld created by the 
neighboring chains. The two order parameters A and h 
have to be found in a self-consistent way: the dimeriza- 
tion amplitude A is proportional to the alternating part 
of ((S„- Sn-fi)), while h is proportional to the alternat- 
ing part of {{S^)), where ((...)) denotes the thermal and 
disorder average. 

Unfortunately, analytical expressions for these disor- 
der averages within the model Eq. (|l|) arc not known and 
their numerical calculation at present cannot be done ef- 
ficiently enough to solve the self-consistency equations. 
This motivates us to simplify the model, by considering 
only the XY part of the spin-spin interaction, i.e., we 
substitute S„-S„+i by S^S^_^_i + -S'^-S'^+i. In the absence 
of disorder the XY model is known to provide a reason- 
able description of both the AF and the SP state ||l2|,|l3|. 
Moreover, as follows from numerical and analytical stud- 
ies, in the presence of disorder the low-temperature be- 
havior of the chain susceptibility to uniform and alter- 
nating magnetic fields is universal, i.e., independent of 
the anisotropy of the spin exchange M-M . 

By means of the Jordan- Wigner transformation the 
XY model can be mapped on a free-fermion Hamiltonian. 
Then, J„^„+i/2 becomes the fermion hopping amplitude 
between sites n and n-t- 1, while h becomes the amplitude 
of an alternating on-site potential. The coexistence of SP 
and AF order in the XY model corresponds to the coex- 
istence of a Peierls dimerization and an on-site charge 
density wave in the spinless fermion model. 



In the weak-coupling and weak-disorder limit, i.e., for 
A,h,SJ <C Jo, we now introduce a continuum descrip- 
tion of the chain (cf . Rcf. [^ ) . The Hamiltonian then 
becomes: 



Hf 



dx^\x) 






a2h 



tp{x), 
(3) 



where aa (a = 1,2,3) are the Pauli matrices. The first 
term in the Hamiltonian Hf describes the free motion 
of the fermions with the Fermi velocity vp = aJo, while 
the second and the third terms describe the backward 
scattering caused by the dimerization, the disorder, and 
the staggered magnetic field. The disorder ?]{x) is related 
to the disorder in the spin-exchange constants by 



r](2na) 



1 



{5J2 



(5J. 



2n,2n+l) , 



(4) 



where a denotes the lattice constant in the chain direc- 
tion. We will assume white noise disorder with a corre- 
lator: 



{ri{x)-q{y)) = A5{x - y), 



(5) 



which corresponds to the statistical independence of the 
variations of the exchange couplings on different bonds 
in the discrete model Eq.(|I]). 

In the absence of a magnetic field {h = 0), the 
disorder-averaged density of single-fermion states p{e) of 
the Hamiltonian Eq.(P) was found analytically in Ref. 
|l8| . The density p{e) is a symmetric function of the 
energy e. Its form depends crucially on the parameter 
g = A/{vpA). For g < 2 the density of states has a pseu- 
dogap (a Peierls gap filled by disorder-induced states), 
while for 5 > 2 (strong disorder) the pseudogap dis- 
appears and p{e) diverges at e = {p{e) ex |e|9 at 

\e\ « A). 

A nonzero alternating magnetic field mixes the h ~ Q 
eigenstates with opposite energies and transforms the 
pair of eigenstates with energies ±e into a pair of eigen- 
states with energies ±-\/e^+7?. Therefore, the disorder- 
averaged r2-potcntial (fi/ = — T(ln2/), S/ being the 
partition function of the grand-canonical ensemble of 
fermions with zero chemical potential) is given by 



2 f^ 
Vlf = --l dep{e)\n 



2 cosh 



(3^/?^+h? 



(6) 



where W is the energy cut-off. The two order parameters, 
A and h, satisfy the self-consistency equations: 



A 



-AA((ai)) 
-A/.((a2)) = 



-Aa 



9A"' 
dh ' 



(7) 
(8) 



where Aa and \h are the corresponding coupling con- 
stants describing the interchain interactions. 

In the absence of disorder, Vif depends on A and h 
only through the combination VA^ H- K^ . As a result, the 
two self-consistency equations acquire the same (BCS) 
form; as, however, they have different coupling constants, 
they cannot be satisfied simultaneously, unless Aa = A?i. 
Thus, in agreement with previous studies [^, we find 
that in the absence of disorder the AF and SP phases can- 
not coexist and the phase with the larger coupling con- 
stant is realized. A competition between these two phases 
always exists in spin chain materials and some special 
conditions, such as a strong spin-phonon couplin g |l3| or 
a significant next-nearest- neighbor interaction |20[| , are 
necessary for the SP state to win. This explains why the 
number of SP materials is small. 

Disorder in the spin-exchange constants suppresses the 
dimerized state by filling the SP gap with single-fermion 
states and thus reducing the energy gain due to dimeriza- 
tion. At the same time, these disorder-induced states en- 
hance the antiferromagnetic susceptibility of the chains: 
The effect of an alternating magnetic field is strongest for 
the fermionic states with |e| < /i, as the occupied state 
with energy — jej is pushed down to —^Je^ -|- K^. The 
higher is the density of states near £ = 0, the more en- 
ergy is gained when AF order appears. Within the mean 
field approximation, this enhancement of the chain mag- 
netic susceptibility due to disorder results in an increase 
of the Neel temperature. 
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FIG. 1. The phase diagram of the disordered SP system de- 
scribed by Eqs.teh,{[iT), and (pi) for A/^ > A^. The dimensionless 
disorder strength A/(i)i?Ao) is proportional to the concentration 
of dopands x (see discussion in the text). The temperature is mea- 
sured in the units of the SP transition temperature at zero disorder. 



From the above we conclude that for A;i > Aa the 
SP state is less favorable than the AF state at all val- 
ues of the disorder strength A. If, on the other hand. 



A A > A;i, a much richer phase diagram arises, as is ob- 
served in Fig. nl This diagram was obtained by numer- 
ically solving Eqs. (|^) and (^ for A a and A/j such that 
Tn{0)/Tsp{0) = 1/4, where T^piA) is the SP transition 
temperature at ft. = and T^{A) is the Neel temper- 
ature at A = 0. Four phases appear: SP, AF, mixed 
SP-I-AF and disordered, separated by second order tran- 
sition lines. At low temperature and weak disorder the 
system is in the SP state. The SP temperature Tsp{A) 
decreases almost linearly with the disorder strength. In 
particular, it can be shown that at small A 



Tsp{A) = TspiO) [1-C 



A 



vpAo 



C = -«1.39, (9) 



where Ag is the value of A for T,A — 0, and 7 = 1.78 . . ., 
is the exponential of Euler's constant. 

Above some critical disorder strength A^f, the system 
undergoes at T/v(A) < Tsp{A) a second (Neel) transition 
into the mixed state, in which the SP and AF orders 
coexist. This coexistence region becomes narrower when 
Xh approaches Aa- Tpf rapidly increases with the disorder 
strength until at some A ^ A^. it becomes equal to the 
SP transition temperature: Tsp{A^:) = Tn{A^:) = T,. 
Above A, only AF long-range order exists and the Neel 
temperature continues to grow slowly with the disorder 
strength. 
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FIG. 2. Detail of the phase diagram Fig. hi. The vertical line 
A/{vpAo) = 0.52 passes through three phase-transition points: the 
SP transition temperature Tsp, the Neel temperature Tjv, and the 
temperature of re-entrance into the undimerized state Tr. Dot- 
ted lines show TSp for A > A* (the SP transition temperature 
calculated at /i = 0) and T^ for A < At (the Neel temperature 
calculated at A = 0). 

The surprising feature of our phase diagram is the fact 
that the disorder strength Asp at which the dimerization 
disappears at zero temperature, is smaller than A*. This 
implies that for Asp < A < A* the system experiences 
three consecutive transitions as the temperature goes 
down (see FigJ^): first the SP transition, next the Neel 
transition, and then the "anti-spin-Peierls" transition, at 



which the SP order disappears. The re-entrance into the 
undimerized state occurs, because the rapid growth of 
the AF order parameter below the Neel temperature sup- 
presses the SP state. 

The existence of a re-entrance transition can be further 
elucidated by considering the Ginzburg-Landau expan- 
sion of il ~ rif + ^ h jT— near the multicritical point 

(r„A): 



fl^aA(T- Tlp{A)) A^ 
+ ah{T^T^{A))h^ + 



&A 

2 ' 



bh ,_4 



■:A^h\ 



(10) 



In Eq.(nG) the coefhcients aA,a/i > 0. Furthermore, 
the stability of the system described by Eq. dlG) , requires 
bA,bh, and D = 6a&/i — c^ to be positive. 

In the presence of a dimerization (at A < A^,) the Neel 
temperature becomes: 



Tjv(A) = r^(A) - — A2. 
ah 



(11) 



As the dimerization suppresses the AF state, c > 0. Sim- 
ilarly, one can find a temperature Tp{A)^ at which A 
becomes zero at nonzero h: 



Tn{A)^np{A)^—h^ 

OiA 



(12) 



which is the temperature of the re-entrance into the 
undimerized state. 

To obtain the dependence of Tm{A) and Tp{A) on the 
disorder strength, we find A and h that minimize fJ and 
substitute them into Eqs. (pi]) and dl2). The result is: 



Tn{A)^T, + {A~A,) 



r dT^j, 


-6Aa.#" 
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(13) 
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Tr{A) ^T, + {A~A, 



From Eq.(|l3|) and the fact that T^{A) increases with 
A, while Tgp{A) decreases, it is easy to find that for 



c> bh 



aA 
ah 



(14) 



both Tn{A) and Tp{A) increase linearly with disorder 
at A < At and Tp{A) < TAr(A). Therefore, inequality 
(O) is the condition for the existence of the re-entrance 
transition. In our model, its validity can be checked an- 
alytically for A/i — > Aa, in which case A^ — > 0, ah — > a a, 
and c — > 2&/j. Our numerical calculations suggest that 
condition (n^) is satisfied for all Xh < \a- 

Next we compare our phase diagram to experimental 
data. At small dopand concentrations x, the observed SP 
transition temperature in Cui-^Zna^GeOa is described by 



Tspix) = Tsp{0) (1 - ax) , 

where a '--^ 14 [g|. To compare this to our result Eq.(^, 
we have to relate the disorder strength A to the dopand 
concentration x. This can be done by assuming that 
the substitution of Cu by Zn changes the spin exchange 
constant by an amount ~ Jq. From Eqs. (0) and M) we 
then obtain: 



(15) the dimerization. These results are in agreement with 
the experimental data on doped CuGeOa. Finally, our 
results indicate the possibility of a re-entrance transition 
from the dimerized SP state back into undimerized state. 
This work is supported by the " Stichting voor Funda- 
menteel Onderzoek der Materie (FOM)" . We would like 
to thank Prof. J.-P. Renard for providing us with the 
manuscript of Ref. H prior to its publication. 



A ~ gJqX. 
Equation (S) then reduces to Eq.(p^) with 



•^0 ^ 



(16) 



(17) 



Keeping in mind that in our model doping-induced large 
fluctuations in some (randomly chosen) spin-exchange 
constants are replaced by small fluctuations on all bonds, 
the agreement with experiment is good. Also the even 
stronger suppression (a ^ 50) of the SP phase in Si doped 
CuGeOa ^^ may be understood: Si, substituting Ge, 
is located between two Cu02-chains and thus influences 
two chains simultaneously |g]. 

One not quite satisfactory aspect of our model is that 
for realistic values of parameters the critical doping con- 
centration Xc at which AF order appears, seems to be 
too large as compared with experiment P| . This may be 
understood from the following considerations: Disorder 
enhances the AF susceptibility of spin chains by filling the 
SP gap with low-energy spin excitations. As was shown 
in Ref. |£l|, with highest probability the excited states 
with energy £ ^ A occur for disorder fluctuations 'ri(x) 
that have the form of a kink-antikink pair. For such a 
fluctuation, the order parameter A(a;) = A + ri[x) has re- 
versed sign in a domain of length R = {vp/A) ln(2A/e) 
between the kink and the antikink. The kink and an- 
tikink, being fractionally charged objects pS], each carry 
spin i which together form a weakly bound singlet. A 
low-energy excited state is then obtained by exciting this 
singlet into a triplet. These weakly bound spins do not 
contribute to the dimerization, but they can give rise 
to AF ordering. However, for weak disorder (small x), 
the density of kink-antikink fluctuations is exponentially 
small in our model, implying that a critical dopand con- 
centration Xc is necessary for AF order to appear. It 
can be argued that Xc would be lowered if the model 
would allow for large fluctuations on some randomly cho- 
sen bonds, in which case the kink density ex x [§|,y|. Such 
an extension of the model, however, makes it much harder 
to obtain the complete phase diagram, as it becomes im- 
possible to make the continuum approximation. 

To summarize, we obtained the phase diagram of a 
disordered SP system, described by a mean-field model. 
We showed that disorder results in a strong suppression of 
the SP state and gives rise to AF long range order, which 
in a certain range of the disorder strength coexists with 
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